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Abstract 

The one-dimensional supersymmetric random Hamihonian -ffsusy = + + where 
4>{x) is a Gaussian white noise of zero mean and variance g, presents particular spectral 
and localization properties at low energy : a Dyson singularity in the integrated density 
of states (IDoS) N{E) ~ 1/ln^i? and a delocalization transition related to the behaviour 
of the Lyapunov exponent (inverse localization length) vanishing like ^{E) ~ l/|ln£'| as 
E ^ 0. We study how this picture is affected by breaking supersymmetry with a scalar 
random potential : H = iJgusy + y{x) where V{x) is a Gaussian white noise of variance 
a. In the limit a <^ ff^ , a fraction of states A^(0) ~ g/\\i'{g^/a) migrate to the negative 
spectrum and the Lyapunov exponent reaches a finite value 7(0) ^ g/h\{g'^/a) si E — Q. 
Exponential (Lifshits) tail of the IDoS for E ^ —00 is studied in detail and is shown to 
involve a competition between the two noises 4>{x) and V{x) whatever the larger is. This 
analysis relies on analytic results for N{E) and 7(£') obtained by two different methods : 
a stochastic method and the replica method. The problem of extreme value statistics of 
eigenvalues is also considered (distribution of the n— th excited state energy). The results 
are analyzed in the context of classical diffusion in a random force field in the presence of 
random annihilation/creation local rates. 

1 Introduction 

The study of spectral and locahzation properties of one-dimensional (Id) random Hamiltonians 
has stimulated a huge activity since the pioneering works of Dyson [T], Schmidt [2], Frisch & 
Lloyd [3j, Lifshits [4j and many others (references to the most important works may be found 
in the review article [5] and the books [SIT]). Because the dimension plays an important role in 
localization problems [8], the strictly one-dimensional situation misses some features of higher 
dimension case (like a weak localization regime). On the other hand the one-dimensional case 
allows to make use of powerful nonperturbative methods and study subtle properties which 
are much more difficult to tackle in higher dimensions. The random Schrodinger Hamiltonians 
-^scalar = + ^{x)-, where V{x) is a random function, have been studied in great detail [B] 

and their properties are rather generic under the asumption that V{x) is correlated on a small 
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length scale and / dx {V{x)V{0)) remains finit(0 : exponential tail in the density of stated at 
low energies (Lifshits singularity) [21 [13l [HI [3l [Til [13 E] and decreasing Lyapunov exponent 
(inverse localization length) at high energy [151 [6] 7 oc 1/E for E — > +00. The situation can be 
quite different if the Hamiltonian possesses some symmetry preserved by the introduction of the 
random potential. Such a situation occurs in the case of supersymmetric random Hamiltonian 

Hsnsy = -^ + Hx? + <l>'{x) (1) 

This Hamiltonian has a positive spectrum, a direct consequence of the fact that it can be factor- 
ized in the form -fTsusy = Q^Q with Q = — ^ + (/)(x) and = 3^ + <P{x)- Moreover, it is worth 
pointing out that i^susy = = Q^Q and its supersymmetric partner H- = QQ^ = + ~ 
(f)' are the two components of the square of the Dirac Hamiltonian Tio = (^2^-^ + o'i(t'{x), where 
(Tj are Pauli matrices : = — + 4>'^ + a^cp' . Therefore the Hamiltonian ([T]) arises naturally 
when studying random Dirac Hamiltonians. Besides its own interest for the physics of localiza- 
tion, this model is relevant in several physical contexts like classical diffusion in a random force 
field (Sinai problem) [M [13 [IS [13 [21 [21] (see §0 below), organic conductors [221 [E] or spin 
chains (the spectrum of excitations of an antiferromagnetic spin-chain is linear at small energies 
like in free fermion model ; the precise mapping of AF spin-chain to free fermions can be achieved 
thanks to a Jordan- Wigner transformation) [231 [24l [25l [26l [271 120] ; see the review provided in 
Ref. [18] . The relation to discrete models has been discussed : the supersymmetric Hamiltonian 
is the continuum limit of a discrete tight-binding Hamiltonian with off-diagonal disorder ^19j . It 
is also the continuum limit of a tight-binding Hamiltonian with diagonal disorder at the band 
center |28l [6] (this point has been recently rediscussed in [29]). The supersymmetry is responsi- 
ble for rather particular spectral and localization properties. For the sake of concreteness, let us 
choose for (f){x) a Gaussian white noise of zero mean, {4>{x)) = and {(j){x)(f){x')) = g 6{x — x'). In 
the low energy limit, E <^ g'^, the integrated density of states (IDoS) presents the Dyson singu- 
larity N{E) ~ 2g/ln'^{g'^/E) [30l [211 [HI [H] , similar to the one of the spring chain with random 
masses [Ij or Anderson model with random hoppings [311 (321 123] . The Lyapunov exponent 
vanishes as 'y{E) ~ 2g/ln{g'^/E) [18] (also obtained for discrete models in Refs. [311 [321 [23] ) , 
dicating a delocalization transition. This delocalization transition is suported by studying other 
quantities : (i) Statistical properties of the zero mode wave function [331 [SI [35] indicate long 
range power law correlations (like the Lyapunov exponent analysis, these calculations do not 
account for boundary conditions), (ii) The distribution of the transmission probability through 
a finite slab of length L at zero energy. In particular the average transmission decreases like 
1/\/L [36j, that is slower than the behaviour 1/L for a quasi Id conducting weakly disordered 
wire, (iii) Time delay distribution presents a log-normal distribution at zero energy [361 [37] . (iv) 
The conductivity is found to be finite at -E = [30l [5]. (v) Finally, the study of extreme value 
statistics of energy levels indicates spectral correlations for E ^ [38j . In the high energy limit 
— > 00, the localization properties are quite unusual since the Lyapunov exponent does not 
vanish but reaches a finite value ^{E — > cxd) ~ g/2. This property is due to the singular nature 
of the potential 0^ -|- (p' with cj) a white noise. When the potential is regularized by introducing 
a small but finite correlation length, it has been shown in [391 [T9] that the Lyapunov exponent 
decreases as 7 oc 1/E for largest energies, as for the random Hamiltonian ^^scaiar = + V{x). 

^ Some interesting results have been also obtained in Ref. [^ in a situation where the correlation function 
growths at large distance like {V{x)V{0)) ~ \x\^ with rj > (the case rj = 1 corresponds to a Brownian motion). 

^ The form of the exponential Lifshits tail depends on the details of the distribution of the random potential. 
Note that the spectrum of the Hamiltonian -ffgcaiar ~ + can also presents power-law singularity : for 

a random potential describing a weak concentration of impurities of negative weights, each trapping a localized 
state at energy Eo < 0, the spectrum presents a power law singularity near Eo, with an exponent proportional to 
the concentration of impurities ; such a singularity is called a Halperin singularity [2] [3l 1101 1111 112] ). 
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If the random function possesses a finite mean value {4>{x)) = fj,g, logarithmic singularities 
are converted into power law singularities \28\ ITS] . Extension to more general situations has 
been considered in Ref. ^40j, where spectrum and localization have been studied for the most 
general random Dirac Id Hamiltonian (random mass, random scalar field, random gauge field), 
however such a study still preserves the (particle-hole) symmetry of the Hamiltonian (note that 
the distribution of the local DoS for this model has been investigated in Ref. [41j). 

The aim of the present article is to discuss the effect of the addition of a scalar random 
potential that breaks the supersymmetry : 



a2 

H = -—r + ^{xf + (l)'{x) + V{x) 
ax^ 



We will mostly consider the case when the functions cp and V are two uncorrelated Gaussian 
white noises with variances {4){x)4){x')) = gS{x — x') and {V{x)V{x')) = a 6{x — x'). The 
case with a finite {(f){x)) will be studied in the section H] with the replica method. The case of 
correlated Gaussian white noises (p and V will be discussed in the appendix|A] where it is mapped 
onto the problem of uncorrelated noises. Our purpose is to study how spectral and localization 
properties of -ffsusy are modified when introducing the scalar potential. A first obvious change is 
that the spectrum of H is not restricted to be positive. Natural questions are therefore : what 
is the number of states sent to by the introduction of the potential V{x), how their energies 
are distributed ? How the delocalization at £^ ^ for the Hamiltonian -ffsusy is affected ? 

The paper is organized as follows. After giving a physical motivation for our model right 
hereafter, we study spectral and localization properties of H in sections [2] and [3] respectively. 
Our approach relies on well-established techniques of stochastic differential equations. In section 
m we employ the replica method in order to find other analytical expressions for the IDoS and 
the Lyapunov exponent and consider the more general case of a finite {(j){x)). 



1.1 A motivation : branching random walks in a disordered environment 

Let us first recall the well-known relation between the Fokker-Planck equation (FPE) describing 
classical diffusion in a force field (/){x) and the Schrodinger equation for a potential </>^ + </>'. 
Let us consider the Langevin equation — = 2(p(x{t)) + \/2r]{t), where the Langevin force 
rj{t) is a normalized white noise. This equation is related to the FPE dtP{x; t) = FxP{x; t) 
where the forward generator reads Fx = — 2dx4>{x). The FPE can be transformed into 
the Schrodinger equation —dtip{x;t) = -ffgusyV'C^sO thanks to the nonunitary transformation 
P{x;t) = ^po{x)■^p{x■,t) since 

Mxr^FxMx) = -Hsusy where Vo(^) = e^' ^"'"^("'^ (3) 

Note that the operator transformation Fx Hgusy is isospectral. ipoix) is annihilated by the 
operator Q defined above : Qtpo = 0. For a confining force field, ipo{x) is the normalizable zero 
mode of i^susy and is related to the stationary distribution of the FPE : P{x;t —>■ oo) ~ ipo{x)'^. 

In order to propose the physical interpretation of the last term of ([2]), we start from a dis- 
crete formulation of the problem of diffusion-controlled reaction in a one-dimensional quenched 
random potential landscape Vk- Let us consider non- interacting particles on an infinite one- 
dimensional lattice with lattice spacing a. We label lattices site by A; E Z, corresponding to 
a position ka. We allow the local occupation number for site k to take arbitrary positive 
integer values (bosonic particles). The transition rates between neighbouring sites k and k + 1 
can be obtained from the Arrhenius law 

tfc+i.fc = ^e^'=-^^-+i (4) 
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where Vk is the potential at site k. The prefactor is chosen in order to obtain a well-defined 
continuum limit a — > 0"*". Additionnally we consider the following chemical reactions : we allow 
particle replication A mA, m ^ 2, with a local rate /3m, fc ^ind particle annihilation A ^ 
with a local rate 7^. The reaction rates are supposed to be random quantities. Therefore, the 
model describes branching random walks in a one-dimensional disordered environment, including 
particle annihilation. 

Let us study the particle distribution on the lattice : we denote the occupation of site k. 
Its mean value obeys the following master equation 

= tk,k+l rik+l + tk,k-l nk-l — (tfe+l.fc + tk-l,k) ^k + {Pk — Ik) rik (5) 

where averaging is taken with respect to the random dynamics defined by rates @ (not to 
be confused with averaging (• • •) with respect to the quenched random potential Vk and random 
annihilation/creation rates). We have introduced fik = Yl°^=i^l^m+i,k- For the continuum 
limit, we introduce the density n{x = ka,t) =nk/a. As a —> we develop \nk±i = n{x,t) it 
adxn{x,t) + ^a^ d'^n{x,t) + ■ ■ ■ Moreover, we introduce the force field 4>{x) via Vk — Vk+i = 
a(j){x = ka) + |a^0'(x = ka) + ■ ■ ■ what allows us to develop the transitions rates dH) as 

1 m <A'(X) 1 Ct>[x) cP'ix) C^ixf 
tk,k±l = ^ =F 7, \ n \ > tk±l,k = -9 =t \ \ \ (d) 

a'^ a 2 2 a'^ a 2 2 

We also introduce the notation 'jk — Pk = V{x = ka) for the difference of annihilation rates 
and creation rates {V{x) > corresponds to annihilation and V{x) < to creation). The 
development yields the partial differential equation 

= ^^^-2-^[cPix)nix,t)]-Vix)nix,t) = -H^^x^t) (7) 

for the average particle density, with H-p-p = —Fx + V{x). We will consider the case where 
the random force field 4>{x) and the random annihilation/creation rates V{x) are correlated over 
small scale. For large scale properties of the diffusion, the minimal model corresponds to assume 
that 4>{x) and V{x) are two Gaussian white noises. The mean value {4>{x)) corresponds to the 
average drift of particles and {V{x)) is related to the average rate of particle annihilation at x. 
We will first consider the case {(pix)) = (the case of finite drift will be discussed in section U]). 
A finite average creation rate {V{x)) corresponds to a trivial global shift of the spectrum of H, 
therefore we will set {V{x)) = 0. 

We have introduced a Fokker-Planck-like differential operator Hpp which, as explained above, 
may be related to the Schrodinger operator ^ thanks to the isospectral transformation Q : 
ip{x,t) = ipQ{x)n{x,t). Hence, the spectrum of H is of great interest for the diffusion problem. 
In particular, if we wish to determine the density n{x, t\y, 0) with initial condition n{x, 0|y, 0) = 
5{x — y) we may rewrite in terms of the spectrum {Ea-, \I'a(x)} of H 

n{x, t\y, 0) = ^ *"(^) ^-(2^) e""""* (8) 

^o(y) ^ 

where il'oix) is the zero mode of -ffgusy given above. A first quantity to consider is the average 
occupation at x at time t after release of a particle at y = x at time t = 0. We can use the 
translation invariance of the problem to identify the position average with averaging with respect 
to disorder : 



I r+L/2 f-+oo 
{n{x,t\x,0)) = lim - / dxn{x,t\x,0) = / dEp{E)e-^^ (9) 

'-L/2 J-00 



L— >CXD L 
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where p{E) denotes the density of states of H (we have omitted averaging in the r.h.s thanks to 
self averaging properties of the density of states). This relation shows that low energy properties 
of the quantum Hamiltonian are related to large time asymptotics for the return probability of 
the classical diffusion problem. 



2 Spectral properties 

In this section we recall the phase formalism, the continuous version of the well-known Dyson- 
Schmidt method [H El [7]. A clear presentation can be found in Refs. [151 E]- The basic idea 
relates on the Sturm-Liouville theorem stating that the number of nodes of the one-dimensional 
wavefunction of energy E is equal to the number of normalizable states below E. The starting 
point is to convert the Sturm-Liouville problen]^ into a Cauchy problerr@ and study the statistical 
properties of the solution of this latter problem. The next step consists to separate the solution 
into an oscillating part and an envelope '4j{x]E) = pe{x) siii6e{x). The study of the phase 
6e{x) permits to analyze the spectral properties of the Hamiltonian H since it allows to count the 
number of nodes of the wave function. The damping of the envelope characterizes its localization 
properties. Strictly speaking, ip{x; E) is the wavefunction only if E coincides with an eigenvalue 
'4j{x]En) oc ipn{x), that is when the second boundary condition is satisfied 'iIj{x = L;En) = 0. 



2.1 Ricatti variable 



It is convenient to start by introducing the "Ricatti" variable z = ip' ftp — (p, the Schrodinger 
equation Hip = Eip leads to the stochastic differential equation (SDE) : 

^z{x) = -E- z{xf - 2 z{x) (l){x) + V{x) (Stratonovich) (10) 

dx 

Since the random functions (p and V are understood to be the white noise limits of some phys- 
ical regular noises (correlated over a finite length scale), the SDE must be understood in the 
Stratonovich sense [32] . The relation (jl20p derived in appendix [B] allows to simplify (jlOp in 
order to deal with one noise only 

dz -{E + z^) dx + ^/a + Agz'^ dW{x) (Stratonovich) (11) 

where W{x) is a normalized Wiener process (primitive of a white noise). We define [3{z) = 
\/ a + Agz"^. This Langevin equation is related to a Fokker-Planck equation (FPE) dxT{z; x) = 
FzT{z; x) where Fz = dz{E + z"^) + i[9^/3(z)]^ is the forward generator. This equation admits a 
stationary solution for a constant flow. The current of z through M corresponds to the number 
of divergencies of the Ricatti variable per unit length, therefore to the number of zeros of the 
wave function per unit length. This is precisely the average integrated density of states (IDoS) 
per unit length N(^E). Therefore 

NiE) = {z' + E) T{z) + \fi{z) ^[/3(z)T(z)] (12) 

We recover on this particular case the general Rice formula lim^^oo •^^^'(-z) = N{E). We intro- 
duce the function U (z) = Ag d 



z 




U{z) =z+ J— ^ - 1 arctan \ \ —z] (13) 



^ A spectral problem is formulated as : find the solutions of H-ip{x) = E-ip{x) for some boundary conditions, 
e.g. ip{0) = 4'{L) = 0. On a finite interval, such solutions E„), exist only for discrete values of the energy 

E £ Spec(H) = {£„}. 

Solve HTp{x; E) — E^{x; E) for given intial conditions, e.g. xpiO; E) = and tp'{0; E) = 1. Solutions exist V E. 
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we obtain the distribution : 



Imposing normalization gives an explicit expression of the IDoS 



(14) 



2.2 Phase and envelope 

The phase formalism introduces another set of variables that give a more transparent picture to 
analyze spectrum and localization. 



Positive part of the spectrum : E = 
= iL + ^(a;). Hi; = Ei; with E 



We write -ffgusy = Q^Q with Q 
: k'^ can be cast in the form 



dx 



k--V{x) ) 



We introduce phase 6 and envelope e^ variables : 



X{x) 



e?(^) sin ( 



cos 9{x) 



+ and 

(15) 
(16) 



(17) 
(18) 



with initial conditions ^(0) = and ^(0) = 0. The phase is related to the Ricatti variable by 
z = — ^ = kcotgO. The interest to deal with this couple of variables lies in the basic idea of 
the phase formalism, i.e. the node counting method : the IDoS coincides with the number of 
nodes of the wave function that can be obtained from the evolution of the cumulative phase. 
The Lyapunov exponent (inverse localization length) is defined as the rate of increase of the 



logarithm of envelope. Therefore N{E) = lima;^oo and 'y{E) = lim^^^oo where we have 
omitted average thanks to self-averaging. These expressions give the most simple way to obtain 
spectrum and localization length from a practical point of view (for numerical calculations). 
Phase and envelope obey the differential equations : 



— = k 



dx 
d£ 

dx 



sin^ 9 + (P{x) sin 20 



V{x) 
2k 



sin26'-</>(x) cos 261 



(19) 
(20) 



Negative part of the spectrum : E = — A;^.- If we perform the same manipulations with 
E = —k"^ , we obtain : 



d9_ 

dx 

de 

dx 



V(x) 

k cos 26 H- sin^ ^ + (t'ix) sin 26 



V(x) 

k sin 26 + — sin 26 
2k 



(x) cos 2(9 



(21) 
(22) 



Invariant measure for the phase.- Using (jl20p we can write (for a positive energy) : 



d6 ' = ' A: dx + P{6) dWix) (Stratonovich) (23) 

where i3{6) = ^ ^ sin^ 6 + g sin^ 26. The related FPE reads dxP{6;x) = FgP{6;x) = where 
Fq = —kde + ^[dg/3{6)]'^ is the forward generator. The current of the phase through the interval 
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[0, vr] is the number of zeros of the wave function per unit length N{E). The stationary solution 
for constant current N{E) = [k - ^(3{e)del3{e)]P{e) is : 

PiO)='-m rj^e""'^ (24) 
13(6) It fSie') 

the IDoS is given by normalizing the distribution. 



2.3 From multiplicative to additive noise 



We have obtained the expression of the IDoS, which is given by normalizing the distribution 
(|14p or the distribution (|24|) and is expressed as a double integral. The analysis of the random 
process and of its distribution is however made more simple by converting the SDE for the Ricatti 
variable ()lip or the phase (I23|) . that include multiplicative noises, into a SDE with additive noise. 
For that purpose we perform the following change of variable : 



dcf <y . , 
z = . — smh 09 

V 45 



(25) 



(that maps M to M). The relation with the phase variable is cotg^ 
variable obeys the SDE 



sinh(^. The new 



a 



Ag 



cosh Lp + \ — 1 



\ (T J cosh if 

where we introduced the potential : 



sinh if + 



dx + v^dTF(x) = -U'{tp) dx + ^/A^dW{x) (26) 



(27) 



AgE 



a 



1 I arctan(sinh y?) 



Note that U{ip) = U{z = sinh 99). 




Figure 1: Potential U{ip) for g = I. Left : a = 0.001 with energy E = 10, 5, 0, -5, -10. Right : 
a = 100 with energy E = 100, 50, 0, -50, -100. 

In order to get the IDoS we contruct the stationary solution of the FPE dxV{ip;x) = 
F^V{ip;x) where F^ = d,^U'{ip) + 2gd'^ is the forward generator. The stationary solution for a 
constant current —N{E) (the variable if goes from +00 to —00, therefore currents for the phase 
9 and for if are opposite) reads : 

V{f) = ^ e-^^(^) r df' e^^(^') (28) 
25 J-00 
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which can also be directly obtained from (jl4p or (j24p since dip = — -v/^ff:^^ = V^t^^- 

An alternative way to obtain the IDoS, that will help the discussion and will be used later, 
is to introduce the n-th moment of the "time" x needed by the process ^{x) to reach —oo, 
starting from (/j(0) = if (spatial coordinate x plays the role of the "time" and variable if of 
the position). This problem is a first exit problem |42j . The moments are given by solving the 
equation B^Tn{ip) = —nTn-i{ip) where = —U'{(p)d^ + 2gd^ is the backward Fokker-Panck 
generator. The solution is constructed for absorbing boundary condition at — oo and reflecting 
boundary at +00 : r„(— 00) = and (9^T„(+oo) = 0. We find (see Ref. [32] or appendix of 
Ref. [38]) 

TM = ^ r d^'e^^('^') /^"d/e-i^(^")r„_i(v.") (29) 

T„(+oo) corresponds to the n-th moment of the time needed by random process 99 to cross M, 
therefore the moment of the distance i between two consecutive nodes of the wave function 
ip{x;E). Let us emphasize on this point. We call £j the distance between the two consecutive 
nodes of the wave function : i/;(0) = ^(^1) = ■0(^i + ^2) = • • • = 0. The problem of first exist 
problem is defined as (p{xi) = +00 and ip{xi = —00 with ip{x) finite for x Gjxj, rcj + The 
random variable is ii and = r„(+oo). Note that all distances are i.i.d. due to the fact that 
potential have a vanishing correlation lengthH. 

The IDoS per unit length is the average number of nodes of tp{x] E) per unit length, what 
corresponds to the inverse average distance between two consecutive nodes : 

iV(E)-i = ri(+oo) (30) 

therefore 

N{E)-' = - d(^e^^(^W dv^'e-^^^^) (31) 
25 J— 00 Jip 



that coincides with the normalization of the distribution (j28p . We will extract limiting be- 
haviours of this exact expression by analyzing more precisely the dynamics of the random pro- 
cess ip{x). 

We first remark that the derivative of the potential at the origin is 




U'{0) = \—E (32) 



For E > the potential is monotonous. 

For E < it developes a local minimum able to trap the process during a finite "time" . In 
this latter case the local minimum of the potential is at ip-^- > and the top of the barrier at 
"P- = ■ 

sinh^j-i- = ±^J (33) 
We easily check that U"{ip±) = ±2y/^ for E < 0. 

Important energy scales.- We will identify later the relevant energy scales in the problem. In 
each reagime {g^ <^ a or g"^ ^ a) two energy scales matter : the two largest scales among a/g, 
cr'^^^, Vgo" and g^. 

• For small supersymmetric noise g^ <C cr, the two relevant energy scales are cr^/^ and cr/g. 

• For large supersymmetric noise g^ ^ a, the two energy scales are ^/go^ and g^. 



^ In general distances it are decorrelated if correlation length is smaller than the length over which deterministic 
dynamics drives ^p{x) to 00. 
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2.4 Density of states for positive energies for ^ cr 



The super symmetric Hamiltonian is characterized by a purely positive spectrum (which follows 
from the structure H, 
[18] : 



susy — Q^Q) which presents the famous Dyson singularity at zero energy p8l 
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for E^O 



(34) 



therefore it vanishes at zero energy : 

j^(a-o)^^ ^ ^ Q_ Y^i^^i ig the fraction of states 
that migrate to M~ when a very small white noise V{x) breaking the supersymmetry is added 



to i?susy ? 



0.15 



^ 0.10- 
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Figure 2: Distribution V{ip) (blue curves) for (a) E = -0.01, (b) E = {c) E = 0.1 and 
{d) E = 2 (with a = 10""^ k g = I). These behaviours can be understood from the shapes 
of the potential of the left part of figure [H For E = —yfog = —0.01 : the distribution is 
compared to exp[— | Y^|^cosh((^ — tf+)]. For the other energies, the distribution is compared 
to approximations obtained in the text. 



Band center : |£'| <^ ^ygcr. In the SDE (|26p . the exponential nature of the potential allows a 
decoupling of the deterministic force and the Langevin force (this works for a only). We 

introduce the value for which the two forces are of the same order : |[/'(<I>o)| == 4(7 : 

^>o ^ In (^leVsVo") (35) 

In the interval [$0 5+c>o[) the dynamics of the random process is governed by the deterministic 
force. The process, starting from +oo, reaches ^>o very fast. Then its dynamic is governed by the 
Langevin force in [— ^>0) '^'o]- Upon arrival at — $o it is driven very fast to — oo by the deterministic 
force. This allows to map the problem to the problem of free diffusion d^p ~ \/^dW{x) on the 
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interval [— ^o]) with reflecting boundary condition at +$o ^-nd absorbing boundary condition 
at — $0 similar approximation was used in p8] to study the supersymmetric Hamiltonian at 
finite energy E <C g^). We immediatly conclude that the average "time" is N{0)-^ = Ti{+oo) ~ 



(distance) 
diffusion 



n:(2$o)^- Therefore a fraction of states 
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iV(0) 



(36) 



have migrated to M^. 

Let us analyze the structure of the distribution V{^p), given by (f28l) . in the low energy limit. 
For (f < — $0 we have ^ 2g therefore jg extremely small and the integral over 

ip' is dominated by the close neighbourhood of : 



^, , N(E) N(E) cosh^>o , ^ ^ 



(37) 



This approximation reflects the fact that when deterministic evolution dominates Velocity ((^) = 
^ ~ —U'{lp) the distribution is V{(p) oc 1/| Velocity (99) |. 

In the intermediate interval [— <l>o,'&o]! e'^^'^)/^^ is almost flat and the distribution is linear 
in this interval 



V{ip) ~ 



N{E) 
29 



1 



,[U{-^o)-U{^)]/2g 



V9 + ^>0 for - $0 < <y5 < +^0 (38) 



where the first term is the contribution of the interval ] — 00, — $0] to the integral ()28p . Finally 
we find for the last interval 



N{E) 
25 



2 2 cosh 99 



for $0 < 9? (39) 



It decreases exponentially : V{(p) ~ e '^+*". The curve is plotted on figure EJ Adding times 
spent in the three intervals gives the normalization A^(0)~^ ~ ^ + i ln^(16Y^ 9^/(^) + '■ 



N{E) ~ 



45 



ln"(2V/t^)+2 



for 1^1 < y/ga 



(40) 



Intermediate energies : y/gcr <^ <^ gf^.- In this limit the potential developes a double 
plateaux structure as suggested on figure [H Once again we use the fact that the deterministic 
force depends exponentially on Lp to decouple the effects of the Langevin force and the determin- 
istic force. The equation |f7'((^)| = 4^ possesses now four solutions : = ±$0 defined above 
and = =b<I>£; with 

^ In {Ej^g) (41) 

The Langevin force dominates the evolution in intervals corresponding to plateaux of U{ip), 
of width — — ln(165(^/-E), while the deterministic force governs the evolution on the 
other intervals. Let us follow the evolution of the process '^{x). (i) In the interval [<I>0)C!o[ 
the deterministic force, dtp ~ —U'{(p)dx ~ — cosh ip dx , drives the process from 99 = 00 to 

93 = $0 in a "time" l/{Ag). (ii) In [<I>£;,<I>o] the Langevin force dominates : dtp ~ ^/AgdW{x). 
Given that one is reflected at <I>0) the average "time" required to reach ^e for the first time 
^d'iffusion ~ '^^^^ ~ ^^)^' (^^^) [—^e,^e] the deterministic force dominates d(/9 ~ 

COSll Ifi 



-U'{ip)dx ~ — w —E — T— dx and drives the process from one edge of the interval to the other 
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in a "time" l/(2g). (iv) In [—^q,—^^]^ the Langevin force dominates : the process crosses 
the interval in an average "time" ^(^o ~ ^e)'^- (v) Finally the deterministic force brings the 
process from —^o to — oo in a "time" 1/(4(7). 

The analysis of the distribution (j28p follows the same logic. In the two intervals where motion 
is diffusive (where the process spends most of the time) : 

, N(E) i ^ + ^0 for -^o<^<-'^E 
V{ip) ~ { (42) 

25 y ^-<^E for ^E<^<^o 

(see figure [2]). Normalizing this distribution gives N{E)~^ ~ g ~'~ 2 ~ In^ {IGg"^ / E) 

therefore we recover the usual Dyson singularity (the scalar potential V(x) plays no role) : 



N{E) - 



ln^(1652/^) +2 



for ^ ^E^g^ (43) 



Large energies : E ^ g^.- Finally, for completeness, we give the distribution in the high energy 
limit. In this case the phase distribution is almost flat P{0) ~ l/vr therefore the distribution for 
99 presents the double peak structure : 

^, , 1 sinh 09+ cosh 03 ,,,, 

V{ip) ~ ^ (44) 

vr sinh ip^ + sinh 

where 99+ is defined by (|33p . The two peaks are associated with inflection points of the potential 
U{(p) where the force is minimum (note that ±(^_|_ ~ ^±)- The IDoS is given by the free IDoS 
N{E) ^^VE. 

2.5 Lifshits tail 

In this paragraph we analyze the tail of the IDoS in the region of rarefaction of states, that is 
for E —00. 

For negative energies, the process (p{x) is trapped by the well at 99 = ip+ a very long "time" 
where positions ip± of the extrema of the potential are given by (j33p . The average "time" needed 
to exit the well due to a fluctuation (Langevin force) is given by the Arrhenius formula. The 
height of the potential barrier is given by : 



lp(..)-.,.,),.;5.(M^) (45, 



with 



F{x) = 2{x + l)[arctan(Va; + 1 + \/x) - 7r/4] - 

r |x3/2 + 0(x5/2) for X« 1 

= <! , (46) 

jx-2^+^+0{x~^/^), forx»l 



Assuming ^[U{(p-) — [/((/?+)] » 1 we can expand integrands in ([3T|) . As we can see on figure 
[H in the limit 5'^ <C cj the potential U{ip) is parabolic near its extrema and we can use formula 
(A22) of Ref. [38]. However in the regime g^ ^ cj the parabolic approximation is not correct. 
In this latter case, noting that 



U{(p) ~ const. ± 2V|£;| cosh(v9 - 99±) , (47) 
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we obtain 




for I-EI > max ( cr^/^ 



(48) 

where Kq(z) is the MacDonald function (modified Bessel function of third kind). We can now 
consider two situations, depending on which among the supersymmetric noise (p{x) or the scalar 
noise V{x) dominates. 

Small supersymmetric noise <^ a - In the intermediate range we recover from (I48p the 
Lifshits tail of the Hamiltonian i^scaiar = --£1 + Vi^) [3 El SH] 



N{E) 



vr 



exp 



for « \E\ « a Is 

OCT 



(49) 



The supersymmetric noise does not affect the DoS in this regime. 
For larger values of \E\ the tail takes the form 



N{E) ~ 



vr 



exp 



+ 2. 
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E\ vr 

3 



5" 



for \E\ > a/g 



(50) 



Eventhough the supersymmetric noise is much smaller than V{x), the behaviour at largest values 
of \E\ is due to a competition between (p and V. 

Large supersymmetric noise a - Expanding (P8|) . we see that the IDoS presents the 

limiting behaviours 



and 



N{E) 



N{E) 



2g 



tt\E\ 



ln\gy\E\) 



exp 



exp 



vr 



tt\E\ 



+ 2 



for ./ga < \E\ < g^ 



for 1^1 > 5^ 



(51) 



(52) 



It is interesting to note that the prefactors coincide with the limiting behaviours obtained for 
positive energies : N{E) ~ i^^^gy^) Vd^ E <^ g"^ and N{E) ~ for E » g"^ . 



2.6 Extreme value spectral statistics 

Up to now we have studied spectral properties through the density of states. In this section 
we consider another property of the spectrum : the problem of extreme value statistics for the 
eigenvalues of the Hamiltonian ([2]). Let us formulate the problem : for a given realization of 
the potential, the spectral (Sturm-Liouville) problem Hiplx) = Eip{x) for boundary conditions 
V'(O) = il^{L) = has a discrete set of solutions Spec(i^) = {En} (we assume that label corre- 
sponds to rank the eigenvalues as -Ei < E2 < E^ < ■ ■ ■)■ We ask the question : what is the 
distribution 

WniE) = {6{E - En)) (53) 

of the n— th eigenvalue ? These distributions give a much more precise information on the 
spectrum than the density of states, what is already clear from the relation Yl'n=i ^niE) = 
Lpi{E) where Pl{E) is the average DoS per unit length accounting for the Dirichlet boundary 
conditions at 2; = and x = L (when L — > 00 the sensitivity to the boundary conditions 
disappears : limi_»oo Pl(-E') = ^'(E), where N{E) is the IDoS per unit length of the infinite 
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system studied above). The distribution Wn{E) gives the probability to find the n-th eigenvalue 
at E whereas the DoS Pl{E) tells us the probability to find any eigenvalue at E. 

The study of extreme value statistics in various contexts has attracted a lot of attention. 
Extreme value statistics of uncorrelated and identically distributed variables were classified long 
time ago (Gumbel for an exponentially decreasing distribution, Frechet for a power law & Weibull 
for distribution with bounded support [44', '45]). Extreme value statistics for correlated variables 
is a much more difficult task. A famous example is the Tracy- Widom distribution for eigenvalues 
of Gaussian random matrices [46 1 Wl\ . There has been a renewed interest in such problems in 
the last years (see for example Refs. [lU H9]). 

The question of extreme value statistics of a Id random Hamiltonian was first addressed 
in Ref. ^SOj for the Hamiltonian H = —-^ + Xln '^n^i^ ~ ^n) where positions are uncorrelated 
and uniformly distributed ; weights Vn are positive, uncorrelated and distributed according to a 
Poisson law. The case of the Hamiltonian //scalar = + ^(^) where V{x) is a white noise 

was studied in Ref. (51] where Wi{E) was derived. This result was generalized in Ref. 
where it was shown that the distributions Wn{E) are Gumbel laws when L ^ oo : despite 
eigenvalues En are random variables a priori correlated, extreme value distributions coincide 
with those of uncorrelated variables. Such an absence of spectral correlations is a consequence 
of the strong localization of the wave functions in this regime [52]. Extreme value spectral 
statistics for the super symmetric Hamiltonian ([1]) near the delocalization transition was also 
considered in Ref. [38] ; it was noticed that in this case the distributions Wn{E) do not coincide 
with extreme value statistics for uncorrelated variables, a consequence of spectral correlations 
near the delocalization transition. 

We first consider the limit of strong supersymmetric disorder g"^ ^ a. When a small scalar 
noise is added to the supersymmetric Hamiltonian, we have seen that the spectrum is not 
anymore constrained to be in M"*". A simple way to obtain the typical ground state energy is to 
write that LN{Ei) ~ 1 from which we obtain Ei ~ —y/o^ InL. Since the exponential tail of 
the IDoS is usually associated to strongly localized states, what will be supported by the study 
of localization in section [3l we expect that the distributions (j53p are similar to the one obtained 
for -ffscaiar (Gumbel laws). This is the aim of the following paragraph to show this statement 
explicitely. 

We first assume that the length of the system is sufficiently long so that the support of (j53p 
is in Mr with energies far from the band center \E\ ^ y/g^- In this case, the process ip{x) is 
trapped by the well of the potential U{ip). The "time" £ needed by the process to go from +cxd 
to — oo {i is the distance between two consecutive nodes of the wave function) is dominated by 
the time needed to exit the well. Its moments ()29p are given by = T„(+oo) ~ n\[Ti{+oo)]^ 
what corresponds to a Poisson law. As a consequence it was shown in Ref. [38] that 

Wn{E) ^ Lp{E) MMTl e-^^(^) (54) 
[n — Ij! 

where the IDoS per unit length of the infinite system is given by (|51|52p . The question of which, 
among (j5T|) or ([52|) . is the behaviour to be considered in order to analyze Wn{E) depends on 
where the support of the distribution is. A priori for the longest size L ^ oo we expect that 
Wn{E) has its support for energies below —g"^ whereas for intermediate length L the support is 
between —g^ and —^Jag. This question will be rediscussed more precisely below. 

We see from eqs. (|51l52p that the density of states per unit length is well approximated 
by p{E) — -^^N(E). In a first time we assume that L is sufficiently large so that the support 

of Wn{E) is below —g^. We can use ()52p from which we write 

Wn{E)c.^—^^(^Te-fi^) (55) 
(n-1)! y/ga \it J 
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with 



n 



m:\E\ 2n\/\E\ vr Fa 



E\L 



■\E\ , 2VI-B 



f{E) = --hi\E\ + '—^ - ' +n-J-^+ ''' ' e " * 



vr 



(56) 



(note that we have reintroduced the term \\foJ^ neglected in ([52]) but present in (j50|) ; this 
will be useful to discuss the other limit <^ a). It is convenient to re-scale energy and length 
as 



TT E 



y 



and L 



go- 



(57) 



and furthermore to introduce the quantity e = {g j'K^ g^)^l^ . Hence, in terms of the new variables 
we find j{E) = g{y)+const. for 



g{y) = n 



with C = Le ^'z")^ 



(58) 



The derivative reads g'{y) = n (1 — ^ — (1 — 2^+2^v^). The first paranthesis vanishes 

for a value of y corresponding to energy out of the range defined in ([52]) ; it should not be 
considered as an extremum. The extremum y = y is solution of 



jO-Vye 



1 



(59) 



In the limit L ^ oo we find : 



lnL + ^lnlnL + 2e\/hl^+ (2 - + O 



In In L e 



nL, 



Rtyp 



vr 



In ( iVh^L 



2^3/4 
^3/2^1/4 



+ 



The width of the distribution is independent on the length : 



(60) 



We can easily show that higher derivatives are given by : g^^\y) — n{—\)^ for L — > cxd and 
k > 1. Typical value of the energies (value that maximizes Wn{E)) is 



(61) 



6E^ ^ 

IT y n 



(62) 



Following |38j we can reconstruct g{y) in the neighbourhood of y by using the derivatives. An 
alternative formulation is to expand g{y — -^^X), where X = j^{E — En^) = ^/n{y — y). We 

can neglect all terms vanishing in the limit y ^ 00 (i.e. L 00) since X ~ 1. We have 
-g (y - -^x) = -lln{y - X/^) + y- Xj^ - 2eJy-X/^ 

+ iJy - X/V^e-S+^/v^+^^V^-^/^ (63) 



Using ([59]) . we obtain 



—g\y ^X ~ const. ^ + e 

n \ yra / y^oo 



n 



(64) 
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Figure 3: nth energy-level distribution Wn{E) for 1 ^ n ^ 6 with g = W, a = 1 and L = 10^. 
We see from eqs. (|61|62p that the dependence in n of the typical energy level and the width are 
E^rl^ ~ E^/^ + i^g^ln(n) and 6En 6Ei/y/E, respectively. 



Therefore we have recovered the Gumbel law 




(65) 



The first distributions are plotted on figure [3l 
Let us do several remarks 



The results (I61I62I65P have been obtained using the asymptotic form of the IDoS (j52p . 
Therefore it was assumed from the outset of the calculation that the support of the distri- 
bution (j65p is below —g^. The condition —E^^ ^ g^ can be recast as a condition on the 
lentgh of the system 



{agy/^ 



(66) 



If the length of the system does not satisfy the condition (j66p . the support of Wn{E) is 
shifted above to the interval between —g^ and —^fag. Therefore the above calculation 
should be redone starting from (j51l) . The results are almost similar : the final distribution 
(j65p still holds for the same width (j62p . Only the behaviour of the typical energy changes 
slightly : 



(67) 




with L' = This expression holds when the length of the system is such that 



9 9 



(68) 



typ 



The crossover (for L ~ e^^/s'^/") obviously corresponds to En 

For smaller system sizes L < 1/g we expect the disorder to have a perturbative effect and 
consequently the energy level to be close to the free levels En — (ttti/L)^, n G N*. FigureH] 
summarizes the different regimes for the groundstate energy Ei. 
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The introduction of the scalar noise has rather strong consequences on the distributions 
Wn{E). (A) For a = distributions Wn{E) are broad distributions (in particular E^'^ ~ 
g^Q-gL g^nd {El) ~ g'^e~^3LY departing from Gumbel distributions, a consequence of 
spectral correlations [38j. (B) For o" 7^ the distributions are narrow distributions centered 
on E^"^ ~ — y^galnL and coinciding with Gumbel distributions, an indication of absence 
of spectral correlations. We now characterize the crossover scale of scalar noise separating 
the two situations (A) & (B). Let us reason at fixed g and L (with L ^ and introduce 
an infinitesimal a : we start from the situation (A). If a is increased, the length L fulfills 
the condition ([68|) and the ground state energy is given by ([67|) , provided that at least one 
state is below —^fgd. This last condition reads LN{—^fga) ~ ^-^i^zj^-^ ^ 1- This leads to 
the crossover value 

' ^ ^ (69) 



~ 



separating (A) and (B). Below this value {a < Uc) the scalar noise can be ignored. Another 



simple way to obtain this scale is to write 

by m- 



^ where the typical energy is given 



Small super symmetric noise.- Finally we mention the results for a ^ . If the support 
of Wn{E) is in the interval between —(j/g and —a^l^ the supersymmetric noise does not 
play any role and we recover the results of Ref. [38J obtained for //scalar : the form ([65]) 
holds for 



ptyp 



3a 



InZ" 



2/3 



and bEn — 



a 



2/3 



2^^ 



InL" 



-1/3 



with L" = ^2nn' (this time the width of the distributions vanish in the limit L 
hypothesis made on the position of the support of Wn{E) implies that the length satisfies 



00 



(70) 



The 



-1/3 



< L < no- 



8 

-1/3 p3 A 



(71) 



For longer lengths we recover the behaviours (|61|62p . Now the term of (|60p neglected 
above is large, what adds a term to (j61|) . We can check that at the crossover (L ~ es V'^/s^) 
both dH]) and ([70]) give E^rl"^ a/g. 

An illustration of the different regimes is given on figure [5l 
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Figure 4: Illustration of regimes for the typical groundstate energy Ei with increasing system 
sizes L in the large supersymmetric noise limit g^ ^ a. 



3 Localization 



Up to now we have concentrated ourselves on the spectral properties of the random Hamiltonian, 
however the most stricking property of Hamitonians with random potentials is the localization 
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Figure 5: Illustration of regimes for the typical groundstate energy Ei with increasing system 
sizes L in the small supersymmetric noise limit <^ a. We recall that C = Le~4\/°vV_ 



of their wavefunctions. In a typical situation, for example if we consider the Hamiltonian 
-f^scaiar = where V{x) is random with short range correlations, one should distinguish 

two regions in the spectrum : in the low energy regime, lowest energy states are those trapped 
by deep wells of the potential. The nature of the trapping depends on the statistical properties 
of V{x) (Gaussian white noise, low density of repulsive or attractive impurities, etc). This kind 
of localization is rather natural. It is correlative to a rarefaction of states reflected in the Lifshits 
exponential tail of the IDoS (section 12. 5p . On the other hand, in the high energy range {E ^ 
disorder) the phenomenon of Anderson localization [53] takes place : in a regime where the static 
potential is a priori perturbative, due to interferences between the extremely large number of 
scattering paths, the wave functions decrease exponentially over distances larger than the Fermi 
wavelength, a nonperturbative effect. Whereas in the 3d situation a delocalization (Anderson) 
transition occurs by tuning the strength of the disordered potential [8] , the Id case is particular 
since all states are localized |54j . a statement rigorously proved in Refs. [55^156). The problem of 
Id Anderson localization has been reformulated and reexamined in many works (see for example 
O O [71 [57] )■ As we mentioned in the introduction the random supersymmetric Hamiltonian 
presents particular localization properties since the low energy Dyson singularity of the IDoS 
[U is accompanied by a delocalization transition [I8]. These features are strongly related to 
the (super)symmetry of the Hamiltonian. In this section we will examine how the localization 
picture is modified by breaking the supersymmetry in the Hamiltonian ([2|). 

Information on localization of wave functions can be obtained by considering different vari- 
ables. The most transparent formulation is probably provided by considering the variables 
{9, ^) of the phase formalism. Localization length £ioc is related to the damping rate of the 
envelope of the wave function. Therefore we can define the localization length by analyzing 
the solution of the Cauchy problem : from Eqs. (I17|18p we take as a definition the relation 
1/Aoc = 7 = lim^^oo where 7 is the Lypaunov exponent (note that we can omit the dis- 
order averaging in this definition thanks to self averaging of this process) It is interesting to 
emphasize that this definition of the localization length is extracted from the solutions ip{x; E) of 
the Cauchy problem, and not from the real wave functions '^n{x) (solution of the Sturm-Liouville 
problem). In other terms the Lyapunov exponent gives a good estimate of the localization length 
of ipn{x) if the statistical properties of the envelope of the solution of the Schrodinger equation 
is not affected when imposing the second boundary condition. In the high energy limit where 
processes 9{x) and ^(x) rapidly decorrelate [15] this is not a problem, however it is not obvious 
that this holds in any situation (in particular for the supersymmetric Hamiltonian -fTsusyj the 

® This picture suggests that the wave function behaves roughly as ipi^) ~ 6*^^ x (oscillations), however one 
should keep in mind that such a simple picture is dangerous since it forgets the important fact that the argument 
of the exponential, ^{x), presents large fluctuations increasing like ^/x (fluctuations of £^{x)/x vanish for x — » 00, 
but not those of £,{x)). The envelope of the wave function is an exponential of a drifted Brownian motion, what 
can have important consequences [58] ; neglecting this important feature can lead to wrong conclusion, like in 
Ref. [59]. 
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Lyapunov exponent does not seem to give a fully satisfactory information as pointed out in the 
conclusion of Ref . [38] ) . 

Since the analysis provided in the previous sections is based on the study of the dynamics of 
z{x) or if{x) = argsinh(y4^/cj 2;(x)), we will extract the localization length from the statistical 
properties of these stochastic processes. We will derive several formulae and use the most adapted 
in the various regimes. Let us recall that the simplest expression of the Lyapunov exponent is 
given by the average of the Ricatti variable p] : 

As in the previous sections we consider here the case (</>) = 0. Together with the expression 
of the stationary distribution T{z), this immediatly gives the Lyapunov exponent. Note that 
since T{z) ~ N{E)/z'^ for |z| — > oo (Rice formula), the expression must be understood as 
7 = J-j^dz z ^) in order to deal with a well defined integral. We can also avoid this 

problem by deriving other formulae, what we do now. 

Positive part of the spectrum : E = We rewrite the two SDE (|19|20p for phase and 

envelope as : 



d(9 = kdx- ^ sin^ 9 dWi (x) + ^ sin 29 dW2 (x) (Stratonovich) (73) 



dC = ^sin20dVFi(x)-^cos20dTy2(2;) (Stratonovich) (74) 

where Wi{x) and W2{x) are two normalized independent Wiener processes. Since the Lyapunov 
exponent is related to {S,{x)) we connect these Stratonovich-SDE to some Ito-SDE and use 
the fact that with this latter prescription, random process and noise are decorrelated at equal 
"time" : 



de = {k+^sm'^9sm29 + ^smA9^dx-^sm^9dWi + ^sm2edW2 (Ito) (75) 



d^ = (^-^sm^9cos29 + gsm^29^dx + ^sm29dWi-^cos29dW2 (Ito) (76) 
We immediatly obtain the following expression : 

7 = ^ = - ^ ^^^^^ ^ + g{sm^ 29) (77) 

where averaging is realized with the stationary distribution. This relation is similar in spirit to 
the one derived in Ref. [60] for //scalar- This equation, with the distribution ()24p . gives another 
explicit expression for the Lyapunov exponent. The Lyapunov exponent can also be expressed 
in term of the distribution (|Ti 



' 2\ (E + z'^Y 

or the distribution (1281) 



(78) 



. = %(WM> with -^M^ ^^^f^-^" im 

Note that the expressions (|77l78l79p are valid for > and are note appropriate to study the 
limit E ^ Q : for example the equation with the Ricatti variable would take the absurd form "7 = 
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'f (^)" (absurd since T{z) is regular at z = 0). The origin of the problem can be understood 



from (1121) that shows that in the limit E 



0, the two terms T{z) = g@ - X^r+m^ 
cannot be considered separately. A more detailed discussion is given in appendix [CJ 



Si^±\R(z)T{z)] 



Band center.- In this regime, due to the previous remark, we start from -f = {z) = {sinh cp) . 
Using the fact that the approximate expression of the distribution is symmetric for \ip\ > ^q, we 
write {sinh If) ~ J^^'^ dipV{(p) sinh99. We obtain 



7(0) 



45 



ln(1653 /a) 



for 1^1 < ./ga 



(80) 



Note however that the multiplicative factor 4 is directly related to our definition of $o separating 
regions where deterministic force and Langevin force dominates : |C/'($o)| = 4(7. Therefore, in 
this derivation, the factor 4 is arbitrary. However the replica method of section U] will predict 
the same prefactor. We won't have the same problem for the other regime since we will use the 
formula ([79]) instead of (172]) . 

This result shows that even a tiny a ^ scalar noise is sufficient to lift the delocalization 
transition of the supersymmetric Hamiltonian. 

Intermediate energies : y/gcr <^ E <^ The function T a{^) presents two symmetric peaks 
centered on ~ ±i ln(4^). Note that \ ln(4A) = \ ln(16E'5'/a) = \{^q + ^e)- We remark that, 
for A = ^> 1, we have d(pT a{^) — 1 and dtp ipT a{'p) — ^ ln(4A) (these equalities 
are already excellent for A = 0.5). Therefore, using the approximate form of the distribution 
derived above, we can write 



2<7 



($0 - ^e) 



/O poo 
dip {ip + $o)Ta(<^) + / dip{ip- $£)Ta((^) 



25 



^>0 - 



Therefore we recover the result obtained for the supersymmetric Hamiltonian alone [18 



25 



ln(1652/^) 



for <^ E <^ 



(81) 



(82) 



High energy limit.- In the high energy limit the distribution of the phase 9 is almost flat, therefore 
using ([771) 

-f{E +00) ~ ^ + I = 73calar + Tsusy (83) 



~ ^ are the high energy Lyapunov exponents for //scalar 



+ 



where 7scaiar ^ ^ and 7susy - 2 
V{x) and Hsnsy, respectively. 

For E ^ 00 the localization length saturates to iioc — 2/g. The high energy wave functions 
present rapid oscillations over a scale 1/k exponentially damped on a larger scale 2/g. 

Negative part of the spectrum : E = —k"^.- As we have seen above, compare to the SDE for 
E = +k'^, the SDE for the variable for E = -k"^ receives an additional term A; sin 20, therefore 



a 



7 = A;(sin26l) - --^ (sin^ 6* cos 26*) + ^(sin^ 26) 



(84) 



In the limit E ^ —00 the phase is trapped at ~ 7r/4 (this is related to trapping of (p by the 
local minimum of potential U{ip) at (p+), therefore 

7(£; ^ -00) ~ + 5 (85) 

This increase of the Lyapunov exponent reflects that the low energy wave functions are sharply 
peaked around deep wells of the potential. 
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4 Replica method 



In this section we derive analytic expressions for the IDoS and the Lyapunov exponent by using 
the rephca method. The computation consists of a shght variant of the method used in Ref. |18j . 
which leads to hypergeometric functions, generalizing the Bessel and Airy functions appearing 
in the pure super symmetric and pure scalar potential problem respectively [61] . Therefore, we 
only sketch the main lines and refer to the Ref. [18] for details. 

We consider the Hamiltonian ^ with V{x) and (/'(x) two uncorrelated Gaussian white noises, 
in the more general case where {4>{x)) is finite : V{x) = ^/ar]{x) and = f^9 ~^ \/^'ni^) ivi^) 
and ry(x) with ^ > are two uncorrelated normalized Gaussian white noises of zero means). 
As mentioned above, the problem of 5-correlations between the noises may be mapped on the 
uncorrelated case (see appendix |Aj) . The spectral properties of H are encoded in the Green's 
function G{x,y;E) given by the matrix element 

a " 

in position space. Here denotes the eigenfunction associated to the energy level Ea, and the 
sum runs over all states a. According to Thouless' formula, average with respect to disorder 
{G{x, x; E)) of the Green's function at equal points yields the derivative of Lyapunov exponent 
as a function of E [62] . Analytic continuation E ^ E — iO+ allows to write 

{G{x,x;E-iO+)) = ^'{E)+mp{E) (87) 

where p{E) is the density of states per unit length. 

4.1 The n-replica Hamiltonian 

We shall make use of the replica trick in order to compute the averaged equal-point Green's 
function (187p (see for example [43]). To this end, we introduce an auxiliary n-component field 
X = {x^, . . . , x") and rewrite {G{x, x; E)) in terms of a Gaussian path integral with respect to x'- 



{{x\{H - E)-'\x)) ^ 



1 /■+^/^ 

- / dx {{x\{H - EY^\x)) 

^ J-L/2 

^lim — / dx Vxxix) (^^P(~2 / dyx{y){H - E)x{y) 
2d d r I ( I \ \ 



-L/2 

Notice that the first line makes explicitly use of translation invariance after average with respect 
to disorder. The limit n ^ eliminates the residual determinant from path integration with 
respect to x- We thus are interested in the n-replica partition function 



= y^x/exp f-^y^ ^ dxx{x){H - E)x{x)\\ = J Vxex.p I- J ^ ^ dxC{x,x) 



\ \ ^ J-L/2 J I \ ■^-L/'i 

where the average over disorder has lead to the Lagrangian 
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As the formula suggests, we abbreviate = the scalar product X ' V — ^iX^V^ 

Rewriting C as 



Vijix) 



9 XiXj 
1 



(91) 



shows that the Lagrangian describes the motion of a point particle in an n-dimensional curved 
space with metric r]ij. The potential is given by 



vix) 



o f X^-^(xY + >)(0)lndet(l+5X^) 

J- + gx z » z 



(92) 



The contact term 6^"'^ (0) may be eliminated by introducting an auxiliary field S = y^l + gx^ and 



rewriting the functional integration measure as Px^^ d{T, 



gx 



1) ; following [18] this term 



will not be considered in the sequel. We recognize an u-model with symmetry group 0(n, 1). In 
one spatial dimension we may transform it to a quantum-mechanical problem in n-dimensions 
where x plays the role of time. Hence we must identify a proper Hamiltonian 7i related to C 
and study its spectrum. TC acts on a Hilbert space with inner product 



(CD, VI/) 



d"xVd^^*(x)*(x) 



d"x 



" a/TT 



^*(x)^(x) 



(93) 



gx^ 



and its eigenvalues £u{n) and eigenfunctions ^u{x) are given by the solutions of Ti.^i,{x) = 
£i^{n)'^^{x) with = {"^u, ^^u) < oo. Since the derivation of H is very much like in [18] we 

only state the result. From ()9ip and (I92p we find the n-replica Hamiltonian 



n 



-A + (l-n)x- V- (x- V)2 + 



n 



1 x' 



2 4 1 + x^ 



(94) 



The angular eigenstates are given by the Gegenbauer polynomials C"^'^ ^(cos 9) where £ denotes 
the main angular quantum number. After separation of the angular part, we are left with the 



radial part of the Hamiltonian that depends only on the modulus p = \/x 

1 d 



2. 



1 ff 



n 



2p dp 2 dp 



ngp d ^ l{l + n 



2p^ 4 1+ gp'^ 



E 



a 



■P 
(95) 



4.2 The groundstate : Lyapunov exponent and IDoS 

In the limit L ^ cxd we expect that the path integral (j89p has a leading term exp[— Lif^G(^)/2] 
where £G{n)/2 corresponds to the ground-state energy of the Hamiltonian TL. Combining (f87l) . 
and (1891) we conclude that 



-f{E) + mN{E) + const. 



2 dZn 



L d' 



'n 



n=0 



dn 



(96) 



n=0 



As above, N{E) denotes the integrated density of states per unit length, and 'j{E) the Lyapunov 
exponent. The constant must be chosen in order to ensure correct asymptotic behaviour E — s- 
iboo (in particular N(E — > — oo) = 0). We shall discuss this problem below. Following the spirit 
of the replica method, we analytically continue £G{n) = n£o + n?£i + • • • and thus identify 



j{E) + mN{E) =£o + const. 



(97) 



21 



We now compute Sq for the Hamiltonian ([95]) . We expect the ground state to be an s-wave state 
with total angular momentum £ = 0. Changing variables to .^^ = 1 + gp^ in (j95p leads to the 
Hamiltonian 



g{e-l) ^ , ng fl d 



• 1 ^ E (y ,^2 



1) + 



1/4 



(98) 



Consequently we must solve the equation Hy"^ = \EQ(n)'^ for the the ground state wave function. 
As for the eigenvalue EQ,{n\ we expand the ground-state wave function into a power series with 
respect to n: ^' = ^'o + + • • • This yields an infinite system of coupled differential equations 
whose first two members are 



T^'n = 



and 



(99) 



Since we seek for a normalizable groundstate wavefunction in the limit n ^ we have to find 
a square-integrable solution of T^'o(0 = 0- Applying the limit ^ ^ 1 in ([99|) . we finally may 
relate ^'o to the eigenvalue 



(100) 



The solution is given in appendix [D] and yields the wave function ^'c 



*o(0 = exp -^,/^ e^'/'U 



1 

+ 2 



E 



a 



a 



4\/^j'^ + ^'TV^ 



(101) 



where U{a,b,z) denotes the second confluent hypergeometric function [63j. Therefore, Eq takes 
the value 




2aU{a + l,fi + 2,i^a/4g^) 



where we have introduced 



2 2 



E 



(102) 



(103) 



/ag 4 V 

The imaginary part can be extracted by using the Wronskian (I14ip of ^'o(0 its complex 
conjugate : 



N{E) 



g_f4£ 

vr V c 



O \ 2 



exp(7rlma) 



|?7(a,/x + l,iv/^7V)|2 



(104) 



We have obtained a compact expression, that can be used more conveniently than the double 
integral (13ip in order to plot the IDoS. 

Eqs. (jl02ll03ll04p provide an exact solutions for the Lyapunov exponent j{E) as well as the 
IDoS N{E) for this model, up to a constant which depends upon a, g and and may be fixed 
by imposing correct asymptotic behaviours, like lim£;^_oo ^^iE) = 0. These results interpolate 
between the known cases of white noise potential and the random supersymmetric Hamiltonian. 

Let us give an example on how to use (|102|103p to study the behaviour at £' = for the 
Sinai case /i = 0. We have 

3 (105) 



1 i 

a = . , ., 

2 SV g-^ 
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tr = 10-' 



-03 -0.2 -0.1 



0.0 01 



(a) 



0.2 0.3 



<T = 10 ■ 
o- = 10~- 



cr =10 



(b) 



Figure 6: IDoS (left) and Lyapunov exponent (right) for g = 1 and /x = for various values of a. 
Delocalization transition at E = for a = (dashed lines) is suppressed even by a tiny scalar 
potential. 



Recall that the confluent hypergeometric function U{a, b, z) behaves like 

i' r(6- l)/^(a)zl-^ b> 1 



U{a,h,z) 



[ (lnz + V'(a))/r(a), 6=1 



(106) 



For small Gaussian noise o" ^ O'^ we tacitely neglect the small imaginary part of a, leading to 
further corrections, and find 



E=o 2\J g 



ln(y^)+W2) 



1 / a 

2y7 



-1 



2g 



. / 1 ro^\ , , ivr , 



(107) 
(108) 



Therefore we obtain the approximate IDoS : 



iV(0) 



In ^g^/a + \n2-il){l/2) +7r2/4 



(109) 



We have recovered by the replica method the behaviour obtained in the sections [2] and [3] 

9 1 /TT, o^ 9 



N{E = 0) 



and 7(^ = 0) 



o--»0 \ii{g'^/a) ' 



(110) 



Note however that the next leading order are different (this is not surprising since the approxi- 
mation scheme of section [2] is quite different). Nevertheless, (|102|103ll04p are less manageable 
for the intermediate regimes singled out in the previous sections. 

Figure [6] illustrates N{E) and ^{E) for the Sinai's case (/i = 0, 5 = 1). Any Gaussian 
noise with cr > lifts the singular bcliaviour N^-^syi^-^') ^ l/(ln£^)^ and. 'ysusy 

{E) ~ l/|ln£;| to 

analyticity in the vicinity oi E = In particular, as shown on figure [U any small a shifts the 
singularity of ^{E) to some minimum at some ii^min > 0. 
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Figure 7: N{E) and 7(-E') for n = 1/4 (top) and fi = 1/2 (bottom). The dashed lines correspond 
to the pure supersymmetric results for £7 = 0. 

The case /x 7^ 0.- It is also interesting to consider the case of finite {(p) = ng. In the absence 
of the scalar noise V{x) (a = 0) the power law Dyson singularity of the IDoS is transformed 
into a power law behaviour N{E) ~ E^^. If a tiny scalar noise is introduced a fraction of states 
migrates to M~ : 



NiE = 0) ^ J- 


"r 




2 

(4?) 


r(^) 



(Ill) 



what we find by straightforward application of (I106P to ()104p . Moreover the feature of smoothing 
singular behaviour extends to < fi < 1/2. For o" = we have the non-analytic behaviour 
7susy(^) ~ ^9 + C± \E\^ with some constants C± for £' > and E < respectively. Again, 
introduction of a shifts this power-law singularity to some minimum of 'y{E) at small positive 
Smim as illustrated on figure [71 In either case, the evaluation of ii^min seems to be difficult. 
However, it would be interesting to find a physical argument for this mechanism. 

5 Conclusion 

In this article we have studied spectral and localization properties of a one-dimensional random 
Hamiltonian H = —■^+(f){x)'^+(j)'{x)+V{x) = Hgusy+yix) which interpolates between the well- 
studied examples of random supersymmetric models -fTsusy and Halperin's model -^scalar- Our 
analysis has pointed out a natural competition between the fluctuations of (j){x) and V{x). We 
have identified the important scales that control this competition for 5^ ^ c or 5^ <C o", which 
are the two largest scales among a/g, o^l'^^ ^Jga and cp' . We have observed that even a small 
additional scalar noise V(x) lifts the singular spectral and localization properties of -fTsusy : the 
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Dyson singularity of the IDoS and the vanishing of the Lyapunov exponent at E = are replaced 
by smooth behaviours : a small additional scalar white noise (a 0) leads to a migration of 
a fraction A^(0) ~ g/\v?{g^ /a) of eigenstates to negative values. It is worth noticing that, V/x, 
the zero energy IDoS (|110lllip for ^ a can be obtained by the substitution E — > ^/ga in the 
known expressions for a = : 

N^"^^\e = {)) ^ n'^''=^\e ^ (112) 

This is a simple consequence of the correct identification of the crossover energy scales. 

Simultaneously to the smoothing of the Dyson singularity, the delocalization transition of 
Hsnsy at = disappears and the Lyapunov exponent takes a finite value 7(0) ~ ^^'^^'^\E ~ 
y^ga) ~ g/ln{g^/a). This logarithmic behaviour shows that, in practice (see figure [6]), even a 
tiny Gaussian noise a kills the singularity of the Lyapunov that becomes almost flat ^{E) ~ g 
for all energies for which density of states is significant. 

IDoS and Lyapunov exponent have also been studied in the other regimes. In particular, 
how the fraction A^(0) ~ g/\v?{g'^ /a) of states are distributed among negative energies has been 
further analyzed ; the precise (Lifshits) exponential tail of the IDoS has been derived in the 
various regimes. It is worth emphasizing that in the lowest part of the spectrum, the tail involves 
a competition between the super symmetric and the scalar noise, N{E — > —00) ~ exp(— 

whatever is the largest scale among g^ (supersymmetric noise) and a (scalar noise). 

The study of spectral properties has been completed by considering the individual distribu- 
tions of eigenenergies (extreme value problem). We have shown that these distributions coincide 
with Gumbel laws, a consequence of the absence of spectral correlations due to the strong lo- 
calization of the wave functions [52], like for the scalar potential alone ; this can be opposed 
to the purely supersymmetric case (cr = 0) for which distributions of eigenenergies are strongly 
modified in the neighbourhood of the delocalization transition [38 j . 

The study of individual distributions of eigenenergies, that includes properly finite size 
(Dirichlet boundary) effects, had allowed us to identify the critical value cJc of the scalar noise 
a below which, for fixed g and L ^ ^/g, the scalar noise can be ignored. We have obtained 
CTc ~ g^e~^^. It is worth noticing that the corresponding value of the E = Lyapunov expo- 
nent (roughly its minimum value) then reads 7(0) ~ g/ ln{g^ /ac) ~ y/g/L. This corresponds to 
a maximum localization length iioc ~ \/L/g <^ L. 

It is not too surprising that the additional white noise modifies spectral and localization 
properties in the vicinity of the band center (around E = 0). However, it is somewhat unexpected 
that, at any value of g (even in the limit g — > O"*"), the noise <^(x) from the supersymmetric 
part controls the spectral properties for E — > —00, what we have seen on the tail N{E) ~ 
exp{—--^\E\) and the distributions of the lowest energy levels. This feature seems counter- 
intuitive since the pure SUSY spectrum is strictly positive so that we would have expected the 
potential V{x) to yield the behaviour N{E) ~ exp(— ^|£'p/^). We attribute this behaviour to 
the singular nature of the supersymmetric potential 0(x)^ + (j)'{x) which is also responsible for 
the saturation of the Lyapunov exponent at high energies j{E) ~ g/2 for E +00. Part of this 
picture will change if supersymmetric noise is replaced by a more regular process with regular 
correlation function of finite width and height (see Ref. [39|). 

Diffusion in random force field with random annihilation/creation rates.- Finally it is interesting 
to come back to the analysis of the results in the context of classical diffusion in random force 
field with random annihilation/creation rates. In order to distinguish more clearly the roles of 
the force field 4>{x) and the annihilation/creation rates V{x), we consider several situations and 
analyze the density of particles {n{x, t\x, 0)) at x at time t, when a particle has been released at 
X initially. Averaging is taken over the random force field and the random annihilation/creation 
rates. 



25 



For g = and <t = : It is useful to recall the obvious fact that in the absence of random 
force field and absorption we have n{x,t\x,0) = ■^/=- 

For g ^ and cr = : Classical diffusion in a random force field (Sinai problem). 
Thanks to ([9]), the spectral Dyson singularity N{E) ~ \/\p? E can be connected to large 
time behaviour [TB] 

(n(x,t|x,0)) ~ -i- (113) 
t^oo In^ t 

much slower than the l/\/t. This behaviour is related to the behaviour x{t) ~ \n^t of 
the typical distance covered by the random walker [18] (see also Ref. [21j where many 
interesting properties of the Sinai problem were studied thanks to the powerful real space 
renormalization group method of Ma & Dasgupta). 

For gr = and cr 7^ : In order to examine the effect of the annihilation/creation 
rates that were chosen to be zero on average, we first switch off the random force field. 
Of course the number of particles is not conserved for cr 7^ 0. In this case the spectral 
Lifshits singularity of the DoS is p{E) ~ ^— exp — ^^^^^ — . The Laplace transform ([9]) 

is dominated by negative energy contributions. A steepest descent estimation shows that 
the averaged number of returning particles diverges with time as : 

1 CT^ 3 

(n(x,i|x,0)) ~ ^=e+48* . (114) 



We emphasize that this increase of the averaged density cannot be compensated by a finite 
mean value of the annihilation rates (y) > that would only add a e~^^^* to this result. 

• For g ^ and <t 7^ : Finally we consider the case of a random force field with random 
annihilation/creation rates. The form taken by the Lifshits singularity p{E) ~ exp— 
leads to the surprising conclusion that the average number of returning particles diverges 
at a finite time tc = n/ ^Jgo : 

(n(x, 0)) = cxD for t tc- (115) 

The two previous points show that this divergence of the average particle density is due to 
the interplay between the random force field and the random annihilation/creation rates. It 
would be an interesting issue to understand precisely the physical origin of this remark. On 
the other hand these last remarks might indicate that the white noise V{x) is probably too 
widely fiuctuating for a reasonnable description of a reallistic random annihilating/creating 
rates. Maybe a more interesting model would be to add a low concentration of such sites. In 
the continuum limit this would correspond to add to the supersymmetric Hamiltonian a scalar 
potential of the form V{x) = '^n'^n^i^ ~ ^n)-, where x„ are random positions with a density 
p and a„ local annihilation/creation rates. The limit of high density p ^ |a„| corresponds to 
the white noise limit studied in the present article. The limit of low density p <C \otn\ might 
be more interesting. This model has been recently studied in the absence of the random force 
field and for absorbing sites (q„ > 0) in Ref. [61], where a penetration length was derived in 
any dimension thanks to renormalization group methods. An interesting question would be to 
understand the effect of the random force field on these known properties. 
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A The case of correlated noises 



It is possible to extent the analysis to correlated noises in the following sense. Suppose that 
V{x) = ^/ar]{x) and 4>{x) = fig + y/ovix) are correlated such that 

{V{x)cl){y)) =T6{x-y) (116) 

We introduce variables C{x) = (j){x)—A and v{x) = 2A{(j){x)— fi g)+V {x) so that the Hamiltonian 
may be rewritten as 



H 



dx^ 



+ C{xY + C{x) + v{x) + 2figA - 



(117) 



The new variables have the correlation function 

{ax)v{y)) = {2gA + T)5{x-y) (118) 
so that the choice A = —T/2g makes them independent. Further characteristics are given by 



{v{x)) = 0, 



{v{x)v{y)) = a5{x- y) 



{C{x)) = ng + —, {C{x)C{y)) = g6ix-y) 
25 



(119) 



where C{x) = — fig — T /2g. Thus, up to a re-definition of energy e = E + fiT + T'^/Ag'^, we 
recover the problem of uncorrelated noises. 

B A useful relation 

Let us consider a random process generated by uncorrelated Wiener processes dWi{t) : 



dx{t) = a{x) dt + bi{x) dWi{t) a{x) dt + V6i(x)6i(x) dW{t) 



(120) 



The equality is valid for Ito and Stratonovich prescriptions. Let us demonstrate this relation. 
Ito's prescription.- Recall that the SDE 

dxi = ai{x) dt + bij{x) dWj{t) (Ito) (121) 

is associated to a FPE dtP = FxP where the Forward Fokker-Planck generator is [42] 

Fa, = -ditti + ^didjbikbjk (122) 

Therefore dxi = a{x) dt + bj{x)dWj is associated to a FPE with generator = —dxa{x) + 
^dxbj{x)bj{x) that is also associated to the SDE dx = a{x) dt + y^bi{x)bi{x) dW{t). Qed. 

Stratonovich's prescription.- The relation between Ito and Stratonovich prescriptions is given 
in Ref. 



dx 



a(x) dt + f3j{x) dWj{t) 



(law) 



dt + PjdWj{t) 



(Stratonovich) (123) 
(Ito) (124) 



dt+ ./p'p'jdWit) 



(Ito) 



(125) 



« + ^/3./3^-^\/M(\/M)' 
a(x)dt + ^ I3j{x)l3j{x) dW{t) 



dt+ y/(5~(5jdW{t) (Stratonovich) (126) 
(Stratonovich) (127) 
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Qed. 

This relation shows that addition law of variances holds not only for additive processes but 
also for multiplicative processes. 



C A remark on the Lyapunov exponent 

In this appendix we clarify some relations between different formulae for the Lyapunov exponent 
given above. 

Let us present the problem with the well-known Halperin model -ffgcaiar = + ^(^)- Here 
V{x) denotes a white-noise potential with average {V{x)) = 0, and {V{x)V{y)) = a5{x — y). 
The widely-used Ricatti mapping allows to relate the spectral statistics for iJgcaiar to passage 
probabilities for a diffusion z{x) whose evolution is governed by the SDE z'{x) = -[E + z{xf] + 
V{x). In particular, the stationary distribution T{z) for z is solution to the differential equation 

^riz) + iz^+E)T{z)=N{E), (128) 

where N{E) denotes the integrated density of states for -^scalar as it can be shown from the node- 
counting theorem. Moreover, the Lyapunov exponent relates to the diffusion via Rice formula 
7 = (z) which, however, must be understood as the principal value 

7= lim / dzzT{z) (129) 



R^ 



in order to avoid difficulties from the asymptotic behaviour T{z) ~ N(E)/z'^ as \z\ -(-oo. For 
E > (fT28]) allows to rewrite 



7= lim [^dzzi^^-—^ ^T'{z)\ (130) 

' R^ooJ_ji \z^+E 2{z^ + E) ^7 ^ ^ 

Clearly, the first term yields 0. Notice that it is crucial to let the integration bounds tend to 
symmetrically, otherwise we would not find a well-defined result. After partial integration of 
the second term, we eventually find an alternative expression for the Lyapunov exponent 

The integration does not require anymore the principal value : it was possible to let the cutoff 
R go to infinity since integrand now vanishes sufficiently fast thanks to the partial integration. 

This relation is particularly useful in order to study the limit E ^ oo since we may immedi- 
ately read of the asymptotic behaviour 7 oc a/E. However the drawback is that (|13ip is rather 
ill-defined for E ^0. Nevertheless, writing 

_ N{E) 2ET{z) aT'{z) 
^^^^"i^^ 2(z2-E)' ^-^"^^^ 

it is not difficult to show that 

/ 2Ez \ a / z'^ + E \ , , 

by partial integration. In order to extract the asymptotics, recall that as E ^ —00 the distri- 
bution T{z) is centered at 2: ~ \J —E. Using this scaling behaviour we recover the asymptotic 
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behaviour 7 oc \/—E. It remains that the E ^ hmit in the two relations (|131|133p seems 
tricky. 

Let us now turn to our model Hamiltonian H = —-^ + 4>{x)'^ + + V{x). Section [2TT] 
provides a detailed account on the Ricatti mapping in this case, in particular the stationary 
distribution T{z) of the variable z{x) was shown to be solution of the differential equation 
N{E) = {z'^ + E + 2gz)T{z) + {a + 4:gz'^)T' {z) /2, see For > we may rewrite 

^>~z^ + E z^ + E 2{z^+E) 
and insert this expression into ()129p what indeed allows to recover ([78]): 

•E+(^^-l] z^\ 

2 \ (^ + z2)2 



7(i? > 0) = ^ ( 7^7— ) (135) 



Conversely, for < we may follow the same strategy as for Halperin's model what yields an 
additional term 

2Ez \ , IE+(^^ + l].^\ 



Again, the advantage of these formulae is that they provide the asymptotic behaviour of 7 as 
E — ^ zboo in a very explicit way. For example, as S — > +00 (|135p shows that 7 oc a/E + 4g 
what is coherent with 7 ~ 7susy + 7scaiar- 

D Solution of the differential equation T^q = 

The differential equation for ^'o is given by 

E a , /i2_i/4 

In the absence of diagonal Gaussian disorder a = we recover the solution given in [18j. For 
cj > we convert the preceding equation into a differential equation for confluent hyper geometric 
functions. Indeed, the ansatz 

^0(6 = exp{\f/2)C"w{z) with z = rj^'^/2 (138) 

where a = 1/2 ± fi, r] = — 2A = ±i^/aJg^ leads to 

w" (z) + (b - z)w' (z) - aw(z) = with 6 = - + a, a = - - — f ^ - ] (139) 

2 2 2r? Vr V/ 

The choice a = 1/2 + n, rj = — 2A = cr/g^ leads to a square integrable solution 



2 .„3« -l) + ^^-lFy^l*oK) = » (137) 



- »p ( -iJZ,.) ^ (^i±i . 1 (^ - 1 7Z) „ , , > /Z:,. ) ,,40) 



In order to see square-integrability, recall that E has a small negative imaginary part —ie. Using 
U{a,b;z) ~ + . . . we find |^'(^)p ~ A second, linearly independent solution is 

z^oo ^-+00 

readily found from complex conjugation of ()140p . The Wronskian which turns out to be useful 
for determination of the integrated density of states may be obtained from known properties of 
U (a, b, z): 

W^(^o(0,^(0) = 2iexp(^Ima) ^'/^-^ exp (^7^?') (141) 
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